Smoothing of Depinning Transitions for Directed Polymers with Quenched Disorder 
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We consider disordered models of pinning of directed polymers on a defect line, including (1 + 1)- 
dimensional interface wetting models, disordered Poland-Scheraga models of DNA denaturation and 
other (1 + d)-dimensional polymers in interaction with columnar defects. We consider also random 
copolymers at a selective interface. These models are known to have a (de)pinning transition at 
some critical line in the phase diagram. In this work we prove that, as soon as disorder is present, 
the transition is at least of second order: the free energy is differentiable at the critical line, and the 
order parameter (contact fraction) vanishes continuously at the transition. On the other hand, it is 
known that the corresponding non-disordered models can have a first order (de)pinning transition, 
with a jump in the order parameter. Our results confirm predictions based on the Harris criterion. 

PACS numbers: 05.70.Fh, 87.15.-v, 02.50.-r 
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Quenched disorder, even in arbitrarily small concentra- 
tion, is expected to modify qualitatively the critical be- 
havior of pure systems in many situations. For instance, 
for Ising spin systems in dimension d < 2 and for systems 
with continuous symmetry and d < 4 it was proven 
via a rigorous version of the Imry-Ma argument |l9l | , that 
randomness in the field conjugated to the order param- 
eter smooths first order phase transitions. An analogous 
result was proven [f| for SOS effective interface models 
in (2 + 1) dimensions. 

In this Letter, we report on a similar phenomenon in 
a very different context, i.e., for (1 + d) -dimensional di- 
rected polymer models interacting with a deject line, and 
for disordered copolymers at selective interfaces. Both 
systems are known to undergo a (de)pinning phase tran- 
sition. Such models have natural applications, e.g., to 
biopoly mers 0, Il2l |26| . to pinning/wetting problems 
|lOl 1 1 11] , to the problem of depinning of flux lines from 
columnar defects in type-II superconducto rs l24l and to 
inhomogeneous surface growth equations |2j|, and at- 
tracted much attention lately, both in the theoretical 
physics and in the mathematical literature. The pure 
(i.e., non-random) models present a variety of critical 
behaviors, ranging from first to infinite order phase tran- 
sitions. On the contrary we prove that, as soon as dis- 
order is present, the transition is always smooth. In a 
way, it is remarkable that one can prove such a general 
result on the nature of the transition, when the knowl- 
edge about the (de)pinning mechanism itself and about 
the location of the critical curve is still quite poor. Our 
result has interesting implications, in particular, on the 
nature of the denaturation transition for inhomogeneous 
Poland-Scheraga models of DNA. 

Pinning /wetting models. We consider general models 
of directed polymers in interaction with a defect line. The 



seemingly abstract setting will be clarified below by some 
physically relevant examples. Polymer configurations are 
sequences S — {S n }n=o,i,... with values in a set E which 
contains a specific point (the origin). We set So = 
0. The free polymer, in absence of interaction with the 
defect line, is described by a homogeneous Markov chain 
on E, with law P. Our only assumption on P is the 
following: let =: r < n < ... be the return times to 
of S (of course, Tj — Ti_i are independent identically 
distributed (IID) random variables). We require 



K(n) = P(n - Ti-i = n) 



(1) 



for some 1 < a < +oo. Logarithmic corrections to the 
power decay are allowed (and actually required for 
a = 1, to make K(-) summable). Note that the first 
return time t% has infinite mean as soon as a < 2. S 
may be transient, i.e. P(ri = oo) > 0. As an example, 
if S is the simple random walk on E = 1 d , then a = 3/2 
for d = 1 and a — d/2 for d > 2. In this case, S is 
transient as soon as d > 3. On the line S = are placed 
quenched IID random charges {w n }n=i,2,... of mean zero 
and variance one. For (3 > and h € M, the Boltzmann 
distribution for the polymer of length N is 
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where H^(S) = En=i(^n - h)l {S„=o} and of course 
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E denoting average with respect to P. Note that the 
polymer-defect interaction takes place only at the con- 
tact points, and that a contact at n with (f3u> n — h) > is 
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energetically favored. On the other hand, polymer con- 
figurations which wander away from the line are much 
more numerous, and therefore entropically favored, with 
respect to those which stay close to it. The main question 
is whether the interaction is enough to pin the polymer 
to the line. The infinite volume free energy of the model, 

F(/3,h)= lim ilogZg,, (4) 

AT— >oc iv 

is self-averaging Moreover, f(/3, h) > 0, as is seen 

restricting the partition function to the configurations 
which do not touch the defect line between 1 and N — 1: 
these paths have zero energy, and their entropy is not 
extensive, in view of JQ. One then defines the pinned 
(or localized) region as 

C={(p,h):F(j3,h)>0} 

and the depinned (delocalized) region as 

V = {(/3,h):F(f3,h)=Q}. 

The denominations pinned/depinned actually correspond 
to the typical polymer behavior. In £ the polymer stays 
close to the defect line and touches it O(N) times before 
the endpoint (various refinements of this statement are 
proved e.g. in |25j for a related model, the copolymer 
introduced below, and more recently in 17] in a more 
general context). On the other hand, in T> the number 
of contacts is at most O(logiV) ^^|. The regions T> and 
C are separated by the critical line h c ((3), so that D = 
{((3,h):h>h c (f3)}. 

The above model has a wide range of applications, and 
a vast literature is dedicated to it. Let us mention two 
particularly interesting examples: 

• (1 + 1) -dimensional wetting of a disordered substrate 
13 El Here, £ = Z+ and a = 3/2. The defect 
line represents a wall with impurities, and S the inter- 
face between two coexisting phases (say, liquid below the 
interface and vapor above), h < means that the un- 
derlying homogeneous substrate repels the liquid phase, 
and vice versa for h > 0. C corresponds then to the dry 
phase (microscopic liquid layer at the wall) and V to the 
wet phase (macroscopic layer) . One of the most debated 
(and still unsettled) issues is whether or not the critical 
line coincides with that of the (exactly solvable) annealed 
model, where disorder is averaged in the partition func- 
tion on the same footing as S. 

• Poland-Scheraga (PS) models of DNA denaturation 
0H3 • In this case £ = Z+, and S n represents the rel- 
ative distance between two DNA strands in correspon- 
dence of the n th base pair: S„ = if the pair is bound, 
S„ > if the bond is broken. Therefore, C (rcsp. V) rep- 
resents the bound (resp. denaturated) phase. Modeling 
S as a simple random walk is known not to be physically 



realistic, and a phenomenological value a > 2 (loop ex- 
ponent), which keeps into account the self-avoidance of 
the two strands, has been proposed j2lj. Therefore, the 
transition is first order in the pure case, cf. © below. Of 
course, real DNA is intrinsically non-homogeneous and 
one resorts very naturally to disordered models like J2J, 
although the IID assumption on u) is very questionable 
in this case. 

Smoothing of the transition. The order parameter as- 
sociated to the transition is the contact fraction, 

f N = JV^E&i (#{1 < n < N : S n = 0}) 

In the pure case ((3 — 0), critical point and critical be- 
havior can be computed explicitly, see e.g. 0,0]. The 
critical point is h c (0) = log(l — P(ti = oo)) < (notice 
that h c (0) < iff S is transient). As for the nature of 
the transition, one distinguishes two cases: it is of first 
order (the contact fraction is discontinuous in the infinite 
volume limit) if X)n>i n K{ n ) < +oo, and of higher order 
if Xm>i n K{ n ) = +oo. In particular, if 5 > then 

f(0, h c (0) - S) ~ const x S for a > 2, (5) 

while 

f(0, h c (0) - 5) - const x tfVCo-i) for 1 < a < 2 (6) 

modulo possible logarithmic corrections. For a — I, the 
transition is of infinite order. 

The main result of this Letter is that, as soon as dis- 
order is present (/3 > 0), the transition is always smooth: 

Theorem 1 For every (3 > there exists < c(/3) < 
+oo such that, for every 1 < a < +oo and 6 > 0, 

F(f3,h c (f3)-5)<ac(f3)5 2 . (7) 

Notice that, since f(/3, h) > 0, (0) is really an estimate on 
the regularity of the transition, an issue debated for ex- 
ample in the context of the disordered PS model 0,0) El- 
ba particular, (0) shows that the order of the transition is 
at least two, i.e., the fraction of bound base pairs vanishes 
continuously approaching h c (j3), in contrast with the con- 
clusions of some numerical studies [l3l Il4| . By convexity, 
self-averaging of F implies self-averaging of the contact 
fraction, whenever 9/ 1 f(/3, h) exists. Theorem^in partic- 
ular excludes the possibility of non-selfaveraging behav- 
ior of the contact fraction at the critical point, which was 
claimed in [lii. Ibl ] . Another interesting consequence of 
Theorem is an upper bound on the number of pinned 
sites in a small window around the the critical point, 
for finite N: indeed, one can show [28J that, if (3 > 
and \h — h c ([3)\ < const x TV -1 / 3 , the probability that 
fN 3> N^ 1 / 3 vanishes for N — > oo. Note that, compar- 
ing © and (7J, our result confirms Harris' criterion 0] 
which, translated into the present context, predicts that 
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disorder is relevant and changes the nature of the transi- 
tion as soon as a > 3/2 (it also predicts that the critical 
behavior does not change if a < 3/2, which is compatible 
with Q.) For previous rigorous work connected to the 
Harris criterion and to critical exponent inequalities for 
random systems, cf. 0. 

As a last remark, note that Theorem ^ is reminiscent 
of the Aizenman-Wehr result Q about smoothing of first 
order phase transitions via quenched disorder in 2d spin 
systems (in particular, the Random Field Ising Model). 
However, the analogy is rather superficial and very dif- 
ferent physical mechanisms are involved in the two cases. 
Indeed Q is based on a comparison between two compet- 
ing effects: on one hand the ordering effect of boundary 
conditions, on the other the effect of random field fluctu- 
ations in the bulk. In our case, instead, boundary con- 
ditions play no role at all (the endpoint Sn is pinned to 
0, cf. (J5J). Our method consists rather in selecting poly- 
mer configurations that visit rare but favorable regions 
with atypical disorder, and in giving Large Deviation Es- 
timates on the number of such regions. This approach 
was partly inspired by where a similar path selec- 
tion method was used to obtain rigorous lower bounds on 
f(/3, h) for the copolymer model. 

A Large Deviations approach. Theorem^is proven in 
full detail in Ref. , under some technical assumptions 
on the law of uj: the result holds in particular if uj n is 
bounded or if it is Gaussian. Here, we present an intu- 
itive argument which clarifies the heart of the method. 
Assume for simplicity a Gaussian distribution for the dis- 
order, uj n ~ Af(0, 1). Let 1 <C I <C N and divide the sys- 
tem into k = N/l blocks B , . . . , Bk-i of length £. For 
a given disorder realization, select the good blocks where 
the sum of the charges is approximately 5£, i.e., let 

l(cu) = I < j < k - 1 : 

[ n=lj+l J 

By elementary large deviations considerations, one real- 
izes that there are typically M typ = (N/£)e~ es I" 1 good 
blocks, two successive good blocks being separated by a 
typical distance dt yp = £e +es / 2 . Next, select all those 
configurations of <S* that touch at the endpoints of the 
good blocks Bj, j 6 1(lo) and that do not touch inside 
the bad blocks Bj,j (/ J(lo) (cf. Fig. 0), and call 
the collection of such configurations. Of course, one ob- 
tains a lower bound on the free energy by restricting the 
partition sum to the selected configurations, i.e., 

^ogZZ > llogE(e<-( S) l {5e ^ } l {Siv=0} ) . (8) 

Thanks to the Markov property of P, the r.h.s. of 
(JHJ) factorizes into a product of terms, one for each 
good block and one for each excursion corresponding 
to a group of adjacent bad blocks (cf. Fig. 0. Note 
that conditioning I independent Gaussian variables to 
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FIG. 1: A typical trajectory in . Here k = 10, £ = 8 and 
I(u) = {3,9}. Note that S n / for n in Bj with j g l(u) 
(bad blocks), except at the boundary with a block Bj with 
j £ I(w). On the other hand, inside Bj, j £ Z(u}) (good 
blocks), the walk moves without constraints. The excursions 
Lo, L\, . . . are typically of length £exp(£8 2 /2). The polymer is 
pinned to zero at steps {0, x,y, z, N}, so that Z^ h ^ factorizes 
into 4 terms. 

have sum 5£ is equivalent, for I large, to shifting the 
mean of each variable from to 5, while keeping their 
variance at 1. Therefore, in each of the good blocks 
the polymer effectively has thermodynamical parameters 
(/?' ',h') = (f3, h — (5/3). Also, note that each of the long 
excursions between two good blocks entails an entropic 
loss \ogK(d t y P ) - -att 2 /2 + 0(\og£), cf. (JTJ. 

Now, take the system at the critical point, h = h c ((3), 
and let N — > oo in (jHJ. By the law of large numbers, 
the free energy contribution of good blocks converges to 
their density, p = £~ x e~ es / 2 , times the average contri- 
bution of each of them, which is £ [f((3, h c (f3) — 5(3) + o(l)] 
for £ large (here and below the error term o(l) denotes 
a non-random quantity that vanishes as £ — > oo). Simi- 
larly, the contribution of excursions converges to p times 
£(~aS 2 /2 + o(l)). In formulas, © implies 

= f(/3, h c ((3)) > e-*£ (f(J3, h c ((3) -50)-?f + o(l)) 

Therefore f(/3, h c (0) —S0) < a5 2 /2 + o(l) for every finite 
I. Since £ is arbitrary, we obtain 0. 

Copolymers at selective interfaces. Consider a poly- 
mer chain close to the interface between two solvents A 
and B, and assume that some of the monomers have a 
larger affinity with A and others with B. If the monomers 
are placed inhomogenously along the chain, the energet- 
ically most favored configurations will stick close to the 
interface. The competition with entropic effects produces 
also in this case a non-trivial (de)localization transition 
at the interface. This model was introduced in the physi- 
cal literature [l2L |2^| and has attracted a lot of attention 
in the mathematical one, cf. e.g. 0,Hj|. The system, 
although physically 3-dimensional, can be reduced to a 
(1 + l)-dimensional one if self-avoidance of the polymer 
is neglected Its Boltzmann distribution can be ex- 
pressed, in analogy with J2J), as 

Pg,(S) = P(5) ^ 1 {SJV=0} , (9) 

Z N,u> 



4 



with the convention sign(O) = +1. Here, the natural 
setting is to take S as a symmetric Markov chain on 
£ = Z, with increments S n — S n _i G { — 1,0, +1}. By 
symmetry, one can take h > 0. The random variables 
((3oj n + h) express the affinity of the n th monomer with 
A, and h is a measure of the asymmetry of the chain 
(if h > there is typically a fraction > 1/2 of monomers 
which prefer A (favorable solvent)). In the literature, the 
only case considered is that of symmetric random walks 
with IID increments S n — S n -\, which implies a = 3/2, 
but in our approach this restriction is not required and 
our analysis covers more general Markov processes. 

Again, one introduces the free energy and, noting that 
in this case f(/3, h) > h/2 (see e.g. 0), one defines local- 
ized and delocalized regions £,T> according to whether 
strict inequality holds or not. Replica method s 1291 and 
real-space renormalization group arguments [23} were 
used to attack the model, and rigorous bounds are known 
for f(/3, h) and for the critical curve h c ((3) separating 
£ and T> 0, 0] . Interestingly, recent numerical simula- 
tions plus probabilistic arguments indicate that none of 
the known bounds is optimal 6], which means that the 
(de) localization mechanism is still poorly understood. 

While the physics of pinning/wetting models and of 
copolymers are rather different, the approach we present 
here is rather robust and works equally well for the two 
problems. Indeed, also for the copolymer model we can 
prove smoothness of the (de)localization transition for all 
(3 > and 1 < a < oc: Theorem[T]still holds, with f(/3, h) 
replaced by f(/3, h) -h/2 [Hj, so that the transition is 
at least second order in view of f(/3, h) — h/2 > 0. Here, 
we give just an idea of how the heuristics above must 
be modified to obtain the result in this case. The main 
point is that © can be rewritten as 



(S) oc P(S)e" 



£„=i(/?« 



hh)A„ 



L {S W =0}> 



(10) 



where A„ = if sign(5„) = +1 and 1 otherwise. In 
this form, the analogy with J2J becomes more evident, 
the role of l/g^o} being played by A„. One can again 
divide the system into blocks and select good ones where 
the sum of the charges is atypically large. However, when 
the selection of trajectories is performed as in JHJ, an 
extra condition has to be met: the selected trajectories, 
S G iS w , must satisfy A„ = for n in a bad block (which 
means sign(SVt) = +1, and not just S n ^ 0). Apart from 
that, the argument is identical as for pinning models. 

Finally, note that for the copolymer the order param- 
eter is no longer the contact fraction /jv, but rather 
f N = A- 1 E JV , W (#{1 < n < N : S n < 0}), i.e., the frac- 
tion of monomers in the (unfavorable) solvent B, which 
vanishes continuously at the transition, in view of our 
result. Again, one can prove finite-size upper bounds of 
order A -1 / 3 for fx, around the critical point |2"sj . 

Conclusions. We have proved that an arbitrarily 
small amount of disorder is enough to smooth the 



(dc)pinning transition in directed (co)polymer models. 
In particular, the transition is always at least of second 
order, even when it is discontinuous in the corresponding 
pure models. Moreover, we have given finite N estimates 
on the order parameter at the critical point. In some 
literature, it is conjectured that the transition is actu- 
ally of order higher than two (possibly infinite) in some 
situations: in particular, for the copolymer and pinning 
models with a = 3/2 [H E3, EH . Our result leaves this 
possibility open. 
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